Quantum stirring as a sensitive probe of ID superfluidity 
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We propose quantum stirring with a laser beam as a probe of superfluid behavior for a strongly 
interacting one-dimensional Bose gas confined to a ring. Within the Luttinger liquid theory frame- 
work, we calculate the fraction of stirred particles per period as a function of the stirring velocity, 
the interaction strength and the coupling between the stirring beam and the bosons. The fraction of 
stirred particles allows to probe superfluidity of the system. We flnd that it crosses over at a critical 
velocity, lower than the sound one, from a characteristic power law at high velocities to a constant 
at low velocities. Some experimental issues on quantum stirring in ring-trapped condensates are 
discussed. 



Progress in the ability to manipulate low-dimensional 
ultracold atomic gases has stimulated the interest in fun- 
damental properties of one-dimensional (ID) Bose liquids 
[11, BH- A Bose-Einstein condensate (BEC) of an atomic 
gas is known to exhibit superfluidity. Experiments have 
confirmed the superfluid behavior by demonstrating a 
critical velocity below which a laser beam could be moved 
through the gas without causing excitations[l,0|, and an 
irrotational flow through the creation of vortices @ and 
vortex lattices in both rotating and nonrotating traps. 
For a Bose-Einstein condensate in a toroidal trap the ob- 
servation of a persistent flow has also been reported [1] . 

Parametric pumping offers another way of inducing 
particle transfer without creating excitations. In pump- 
ing, periodic (ac) perturbations of the system yield a dc 
current. Indeed, this current may be entirely adiabatic 
as long as the external perturbations are slow enough 
such that the system always remains in the instantaneous 
ground state. The number of particles transferred in each 
cycle is then independent of the pumping period T and 
the integral of the current over a period is quantized for 
a clean inflnite periodic system with filled bands 0, [lo| . 
Up to now, spectacular precision of quantization of the 
pumped current has been achieved in experiments with 
nano-electronics devices 

Quantum pumping is intimately connected to quan- 
tum stirring. Quantum stirring is accomplished by the 
cyclical variation of one system parameter, while preserv- 
ing the characteristic of a pump, i.e. the orientation of 
the particle flow is fixed. Quantum stirring has been ex- 
ploited to elucidate the nature of the critical velocity in 
superfiuid liquids 0, @, [l2|- We focus here on stirring a 
one-dimensional (ID) interacting Bose gas with a laser 
beam in the regime where interaction effects are espe- 
cially strong, and we propose the fraction of stirred par- 
ticles as a measure of superfluid behavior. Although for 
a homogeneous ID Bose gas the superfluid fraction, de- 
fined as the response to twisted boundary conditions and 
estimated from ground-state quantities (see e.g. is 
always one independent of the interaction strength, it is 
a relevant question to ask whether the out of equilibrium 



behavior of a ID Bose gas under the effect of an external 
probe is closer to the behavior expected for a superfluid 
(e.g. frictionless flow below a certain velocity threshold) 
or rather to a normal fluid (e.g. flow with drag). The 
study of the stirred fraction gives a measure of the degree 
of superfluidity of the fluid, i.e. a small stirred fraction 
corresponds to superfluid-like and a unity stirred frac- 
tion corresponds to normal-like behavior. This is com- 
plementary to the onset of a drag force as manifestation 
of superfluid behavior 13j, [14 1 . 

We consider N bosons of mass m confined onto a 
ID ring of circumference L, with contact interactions 
v{x — x') — gS{x — x') at zero temperature. The 
long-wavelength behavior of this system at distances 
larger than the cutoff length a = l/po = L/N is de- 
scribed by the Luttinger liquid Hamiltonian in terms of 
the density and phase fluctuation modes of the bosonic 
field [lllii[i3l: 
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the field 9{x) corresponds to the phase of the super- 
fiuid, and we have [(/)(a;), V0(a;')] = iS{x — x'). In the 
case of repulsive contact interaction between bosons, 
the Lutt ing er parameters Vs and K used in ([T]) are 
obtained [1711 by: VgK = ^^^7^, as follows from galilean 
invariance, and ^ = in the weak coupling limit, 

while ^ = 1^ (1 - in the strong coupling limit. 

When the interaction goes to zero K goes to infinity, 
while K ^ \ for infinitely strong hard-core interac- 
tions (Tonks-Girardeau limit), where the problem can be 
solved by mapping onto a gas of noninteracting fermions 
[ist - In this regime 27rpo — > 2fci?, with kp being the 
Fermi wavevector of the corresponding mapped spinless 
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fermions. The long-wavelength properties of ID dipolar 
gases are also described by Eq. ([J) with K < 1 [l9|] . 

We describe next the effect of a barrier moving with 
velocity V through the fluid by introducing the time- 
dependent potential U{x,t) = UoS{x — Vt). In an ex- 
periment this could be realized eg by stirring the gas by 
a blue-detuned laser beam. The Hamiltonian acquires 
an explicitly time-dependent term which couples to the 
density: 



SH{t) = / dxU{x,t)p{x). 



(3) 



Using Eq. ^ for the density and keeping only the lowest, 
most relevant harmonics we may rewrite ([3]) as 

SH{t) = Uo[pa - -V0(W) -f- 2pQ cos(27rpoV"i - 2(f>{Vt))]. 

TT 

(4) 

The term proportional to V0 is analogous to a slowly 
varying chemical potential and can be absorbed in H by 
a redefinition of the field 4>t 4> ^ 4>— {K/vg) dx'U{x'), 
while the last leading term in Eq. ^ represents scatter- 
ing of the bosons off the barrier with momentum close to 
±27rpo- In the Tonks Giradeau^lSj limit it describes the 
backscattering of right-movers into left-movers, i.e. pro- 
cesses with momentum close to ±2fci^. During its motion 
the barrier drags along a part of the bosons. We are 
interested in the stirred fraction Ngtir/N i.e. the frac- 
tion of particles transported per period T = L/V by 
the moving barrier, and related to the particle current 
as A^stir = ^ dtl{t). If the barrier height is infinitely 
large, the fraction of stirred particles per period is quan- 
tized, i.e. Nstir/N = 1, independently of the interaction 
strength. If the barrier height is finite, the stirred frac- 
tion is in general smaller than one and we show that it 
is related to the degree of correlations in the system. We 
analyze perturbatively the regimes of weak and large bar- 
rier for arbitrary interaction strength and treat exactly 
the Tonks-Girardeau regime. 

Weak barrier- In the weak barrier limit we perform 
a perturbative analysis of the current generated by the 
stirring Hamiltonian SH. As customary in Luttinger 
liquid formalism we introduce the particle density of 
right(left) movers related to the fields 9{x) and as 
Pui^L-) ^ ^ ± [V0(a;) =F V(/)(a;)]. The particle current at 
low energy is J{x) ~ V9{x); since it involves the diff'er- 
ence in the number of right and left movers, the term 
proportional to in Eq. (j4|), which does not distin- 
guish between left and right movers, plays no role in gen- 
erating the particle current. On the contrary, the third, 
backscattering term in Eq. ^ , can lead to the generation 
of a current which we define of backscattering, lb- In fact, 
addition of the moving-barrier potential breaks the con- 
tinuous chiral svmmetrv(2l| violating the conservation of 
the axial charge — Nl, where Nj^^i^) = J dxpii(L){x)- 
In the lowest order perturbation theory the backscatter- 
ing current is given by /° = I [iV^ , (5iJ] = -jr\Nu,5H]. 



In our specific case, by using the bosonized expres- 
sion of the density operators and of the stirring Hamil- 
tonian, the resulting backscattering-current operator is 
= i-f{t)n{t) - h.c, where -f{t) = C/oe^^'^PoVt ^nd 
n ~ poe*^'^^^'\ and it is characterized by the backscat- 
tering frequency ujh = 2TrpQV. Linear response theory 
yields the backscattering current to second order in the 
barrier strength Uo as 



dt'{[l',{t),SH{t')])H, 



(5) 



and turns out to be related to the Fourier transform of the 
Green's function of the backscattering operator e*^'^^^*) 
at the characteristic frequency Wf,. In the thermodynamic 
limit N, L —>■ oo with pq = N/L constant and for small 
stirring velocity, the resulting backscattering current is 
given by 
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with r being the Euler Gamma function. The fraction of 
stirred particles is readily obtained from the backscat- 
tering current according to Ngtir/N — Ib/oJb- In the 
Tonks-Girardeau hmit K ^ 1, Eq.® yields Nsta/N oc 
(Uq/Hvs)'^, i.e. the result is independent of the frequency 
ujb and hence adiabatic In the small ujb limit this 
result is in agreement with the exact calculation of the 
fraction of stirred particles, as shown below. Note that 
as the Luttinger liquid theory is an effective low-energy 
model, it describes correctly the system at frequencies 
ujb < l-RVsja, hence the expression ([6]) is valid only if 
V < and cannot treat the supersonic regime. By 
recalling that the power-law dependence in Eq. ([S]) origi- 
nates from the excitation of sound waves in the quasi-one- 
dimensional geometry, we can also determine the smallest 
velocity for which Eq. ([6]) holds in the case of a ring of fi- 
nite length. In this case no excitations are possible below 
the lowest velocity Viow = Vs/N ~ nh/mL corresponding 
to the momentum of the lowest bosonic mode on the ring. 
The value of Viow found agrees with the one obtained by 
using a Gross-Pitaevskii approach for K ^ 1 Thus 
as a main result we find that at the critical velocity Viow 
the fraction of stirred particles crosses from a power-law 
to a constant (adiabatic regime). Note also that the adia- 
batically stirred fraction decreases with decreasing inter- 
action strength as l/r{2K): when K grows, the system 
becomes more superfiuid, hence the interaction with the 
external barrier decreases and Nstii/N — > 0. 

The results obtained above are consistent with a treat- 
ment based on the perturbative renormalization group 
(RG) approach 2j|. In this approach the scaling of the 
potential Uq with frequency uj is obtained from the flow 
equation dUo/dz — {K — \)Uo where dz = duj/cu. As 
a function of K, two regimes are distinguished. When 
K > 1 the barrier is irrelevant: Uq decreases as lu is 
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FIG. 1: (Color online) Left panel: Fraction of stirred particles 
as a function of the stirring velocity (in units of Vs), obtained 
matching Eqs.(|6]l and dTjl at Vy^ = Vt^ = V* through the RG 
considerations (see text). (Black) Continuous line: Nstii/N at 
decreasing barrier strengths Uo/(hvs) from 1. (upper curve) 
to 0.2 (lower curve) for K > 1, fixed as 2; (red) dashed line: 
Nstir/N at increasing tunneling strength to / (pohva) from 0.2 
(upper curve) to 0.6 (lower curve) for K < 1, fixed at 0.5. 
The vertical dashed-dotted lines indicate V* /vs- The vertical 
straight line indicates the critical velocity Viow below which 
the stirring is adiabatic. Right panel: Summary of the RG 
flow for the barrier potential Uo and tunneling strength to 
(left and right edge arrows) and Nstiv/N (arrows in the middle 
of the frame) at varying the interaction K approaching the 
adiabatic regime. The stirred fraction is the analogous for 
neutral particles of the conductance for a ID electron gas 
with a barrier [2^ 



In the presence of tunneling the stirred fraction of parti- 
cles is Nstir/N = 1 — It/uJb, where It/oJt is the fraction 
of tunnelled particles, not stirred. In the hard-core limit 
{K — 1) the stirred current will be again linear in the 
frequency of the stirring. We thus recover the adiabatic 
limit. Under the RG flow, the tunneling becomes relevant 
for interacting bosons with contact repulsion {K > 1) 
therefore, upon decreasing the stirring velocity the effec- 
tive tunneling strength increases, thereby decreasing the 
stirred particle fraction, again shown in the right panel 
of FiglU Perturbation theory breaks down when the ef- 
fective tunneling strength reaches unity and the RG flow 
must be stopped at V^* = Vs{to/pohvs) , then the 
stirred fraction of particles is governed by the previous 
weak barrier limit. The results for the dependence of the 
stirred fraction of particles on the velocity V is shown 
in Fig. [T] The results explicitly show a difference in the 
regime with K > 1 (short-range interactions) and K < 1 
(dipolar interactions). In the latter case the stirred frac- 
tion decreases at increasing velocities, where the system 
tends towards superfluid behavior. Since a dipolar gas is 
characterized by a quasi-crystal order phase at increas- 
ing density yjl], the result can be interpreted as an in- 
efficiency of the stirring in creating an excitation in the 
ordered state. 



decreased from Vs/a down to ujb ~ V/a. For an infinite 
system, Uq and hence Ngtii/N scale to zero as ujb — * 0; for 
a finite system the RG procedure should stopped when 
Vs/ L, i.e. for V ~ Vs/N; this is the regime where we 
find a residual adiabatically stirred fraction, independent 
of V. For K < 1, e.g. in the dipolar gas, Uq and hence 
-^stir/-^ grow under RG, i.e. the barrier is a relevant per- 
turbation. This is shown in the right panel of Fig[TJ Per- 
turbation theory breaks down when NsUr/N ~ 1, i.e. at 
the velocity = Vs{Uo/hvs) ^-'•^ and the RG flow must 
be stopped. The behavior beyond this breakdown point is 
described by an effective weak-link tunneling model 24 1 . 



Weak link limit- The large-barrier limit is equivalent to 
a ring cut at the position of the delta barrier, and we treat 
the residual tunneling between the two ends of the ring 
as a perturbation (2 4] . In this case the bosonized Hamil- 
tonian corresponding to the hopping across the weak-link 
can be obtained by a duality transformations ip ^ 9 
and is given by SH ^ to cos{2TTd{Vt)). Its contribution 
to the particle current (tunneling current It) can be cal- 
culated in the linear response regime and its explicit ex- 
pression for an infinitely long ring is: 



It = 



{pohvsfT{2/K) \v. 



— 2ttpoV. 



(7) 



Non-perturbative analysis- In the Tonks-Giradeau 
limit {K = 1) a time-dependent Fermi-Bose (FB) 
mapping flst is employed to generate exact solu- 
tions of the problem[25| and the current is calculated 
exactly. The time dependent version of the FB 
mapping permits to write the exact many-body wave- 
function of N impenetrable bosons on a ring as 

'^b{xi,. . . ,XM\t) = A{xi, . . . ,XN)'^Fixi, . . . ,XN:,t), 

where A is a unit antisymmetric function 
A{xi,...,xn) = ni<j<fc<Ar sgn(a;fe - Xj), and 
^'^(xi, . . . , xat; t) — C det^j^i'ipi{xj,t) is the wave 
function for any ideal Fermi gas, ^i{xj,t) be- 
ing the solutions of the one-body time-dependent 
Schroedinger equation in the external potential U{x,t). 
Starting from the above many-body wavefunction, 
we evaluate the Tonks- Girardeau particle current 
density in terms of the one-body density matrix 
pi{x,y) = j dx2.-.dxN'^*B{x, - ■ ■ ,XN]t)'^B{y, - ■ ■ ,XN]t) 
as J[x) = —{h/2mi)[drPi{x-\-r/2,x — r/2)\r=o. Although 
pi{x,y) for a TG gas is very different from the one of a 
Fermi gas due to the presence of the mapping function 
A, we find that the latter has no effect on the current, 
which then coincides with the current of an ideal Fermi 
gas. In the adiabatic limit V < nh/mL the particle 
current and the stirred fraction produced by the slow 
variation of the stirring potential can then be evaluated 
by following the adiabatic expansion of Thouless for an 



4 




FIG. 2: Fraction of stirred particles for a Tonks-Girardeau 
gas (i.e. case if = 1) as a function of the barrier strength 
Uo/{hvs), in the adiabatic hmit from Thouless expression 
Eq.(|9]) (solid line) and analytical small-[/o behavior (dot- 
dashed line) as explained in the text. 



ideal Fermi gas[9|, i.e. 



27r mL 



dt 



E 



[{i>j\Tpi) {dxil^i\ipj)+h.c 
^ (8) 

where fij is the fermionic probabiUty occupation func- 
tion of the state In the case of a blue-detuned laser 
field piercing the ring at a position x — and modelled 
by the potential UoS{x), the appropriate orbitals ipi(xj, t) 
are the L-periodic free-particle energy eigenstates satis- 
fying at X = the cusp condition. The complete or- 
thonormal set of even parity V'i^^ and odd parity ipn 
eigenstates are ip^r^\x) = (e'''="^ -|- e~'''"^'^~^'>)/Afn and 
tpn \x) = 'ij^ sm{2mTx / L) , where fe„ are obtained from 

the trascendental equation fc„ tan(fc„L/2) = mllo/fi'^ (for 
Uq oo we have k„ = 7r(2n + 1)/L, in agreement with 



jlS'l) and A/'n = ^/2L\T+ sin(fc„L)/fc„i)], with n running 
from 1 to CX3. The A'^-fermion ground state is obtained 
by inserting the lowest N orbitals into the determinant 
above (Fermi sea), and using the exact orbitals V""^^ as 
instantaneous ground state we obtain from ([8|) 



N. 



= 64^(/,-/,)- 



sin2(fcjL/2) (fcjL)247r2^2 



sm{kjL)/kjL [{kjL)^ - Air^e^]^ 

(9) 

For a weak barrier by using the small-?7o expression for 
kn we obtain Ng/N ~ 0.32{Uo/hvs)^, which scales as the 
K = 1 limit of the backscattered current in Eq.® be- 
cause Vs — kkp/m for K = 1 (2^ . For an infinitely strong 
barrier using the J7o — > oo limit of fc„ it's straightfor- 
ward to verify that the particle transport is quantized 
i.e. all the particles are dragged by the barrier and 
./Vstir/-/V — 1. This is shown in FigllJ where the stirred 
faction of particles is plotted as a function of the barrier 
strength. 

Experimental issues on condensates in closed loop 



waveguide -A possible way of achieving experimentally 
an annular condensate with strong transverse confine- 
ment is to use a magnetic toroidal trap, as reported in 
[3, 111] ■ Experimentally the stirring of hydrodynamic 
flow in a BEC by a blue-detuned laser beam, has been 
analyzed by calorimetric method^] and phase contrast 
imaging^. The onset of a drag force has been shown 
by the asymmetry in the density profile, defined as the 
difference between the peak column density in front and 
behind the laser beam, as a function of the stirring ve- 
locity above a critical velocity. The space integral of the 
density asymmetry is analogous to the fraction of stirred 
particles calculated above. Recently, the persistent flow 
of Bose-condensed atoms in a toroidal trap has also been 
observed 0. A variant to such experiment by the addi- 
tion of a cyclic moving plug beam could be a valuable 
realization of the present proposal. 

In conclusion, superfluid flow in a ring geometry raises 
interesting new possibilities. With the use of a moving 
barrier acting as a quantum stirrer, the analog of quan- 
tization of particle transport for electron systems could 
be realized for a gas of atoms as an alternative probe of 
superfluidity. 
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